How pure is the tail of gravitational collapse ? 
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Abstract 

Waves propagating in a curved spacetime develop tails. In particular, it is well established 
that the late-time dynamics of gravitational collapse is dominated by a power-law decaying tail of 
the form Mt~( 2l+3 \ where M is the black-hole mass. It should be emphasized, however, that in 
a typical evolution scenario there is a considerable time window in which the signal is no longer 
dominated by the black-hole quasinormal modes, but the leading order power-law tail has not 
yet taken over. Higher-order terms may have a considerable contribution to the signal at these 
intermediate times. It is therefore of interest to analyze these higher-order corrections to the 
leading-order power-law behavior. We show that the higher-order contamination terms die off 
at late times as M 2 t~ 4 ln(i/M) for spherical perturbations, and as M 2 t~( 2l+4: ^ In 2 (t/M) for non- 
spherical (I ^ 0) perturbations. These results imply that the leading-order power-law tail becomes 
"pure" (namely, with less than 1% contamination) only at extremely late times of the order of 
10 4 M. 
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One of the most remarkable features of wave dynamics in curved spacetimes is the de- 
velopment of "tails". Gravitational waves (and also other fields) propagate not only along 
light cones, but also spread inside them. This implies that at late times waves do not cut off 
sharply but rather die off in tails. Price [l| was the first to analyze the mechanism by which 
the spacetime outside of a (nearly spherical) collapsing star divests itself of all radiative 
multipole moments, and leaves behind a Schwarzschild black hole. 

It is well-established that the late-time dynamics of massless fields propagating in the 
curved Schwarzschild spacetime is dominated by inverse power-law decaying tails, the power 
indices equal 21 + 3 [lj. Physically, these power-law tails are associated with the backscat- 
tering of waves off the effective curvature potential at asymptotically far regions E ] . 

The analysis of Price has been extended by many authors. Gundlach, Price, and Pullin 
j^l showed that power-law tails are a genuine feature of gravitational collapse- the existence 
of these tails was demonstrated in full non-linear numerical simulations. Furthermore, since 
the late-time tail is a direct consequence of the scattering of the waves at asymptotically 
far regions, it has been suggested that power-law tails would develop independently of the 
existence of an horizon [4|. It is worth mentioning that the phenomena of wave tails 
been analyzed in the non-spherical (rotating) Kerr spacetime [5, 
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Despite the flurry of research in this field, our current understanding of the evolution 
of wave tails in curved spacetimes is somewhat unsatisfactory. The leading order inverse 
power-law tails are well established, but the resultant formulas are only truly useful at very 
late times. In a typical evolution scenario there is a considerable time window in which the 
signal is no longer dominated by the black-hole quasinormal modes, but the leading order 
power-law tail has not yet taken over 19j. In this intermediate asymptotic regime (which, 
as we shall show below, may be of a very long duration) the signal may have a significant 
contribution from higher-order terms. These correction terms actually contaminate the 
leading-order power-law tail. 

The analysis of higher-order corrections to the leading-order power-law tail is of theoret- 
ical as well as practical importance. This is especially crucial for the determination of the 
power index in numerical simulations which are always characterized by finite integration 
times. The precise (pure) power index is expected only at infinitely late times. However, 
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in any realistic (finite time) numerical simulation there would be higher-order terms which 
would contaminate the pure inverse power-law tails. The influence of these higher-order 
terms becomes more significant as the available time of integration decreases. Thus, in 
practice, the limited time of integration introduces an inherent error in the determination 
of the power index. It is therefore of great interest to calculate the higher-order correction 
terms. The ratio of these contamination terms to the leading order term would then pro- 
vide a systematic indication of the error introduced by the finite-time numerical simulation. 
In this work we generalize our previous calculations 26| to include generic non-spherical 
perturbations. 

We consider the evolution of a massless scalar field in the spherically symmetric 
Schwarzschild background (a collapsing star of a fixed black hole). The external gravi- 
tational field is given by 



ds 2 
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Resolving the field into spherical harmonics $ = ^2i lm ^ l m (t, t)Y^ 
^;ab9 ab ^ — for the scalar field in the curved background reads 

Q2 Q2 



dr 2 + r 2 dQ 2 . (1) 
, 4>)/r, the wave equation 



+ V(x) 







(2) 



.dt 2 dx 2 

where the tortoise radial coordinate x is defined by dx = dr{l — 2M/r). The effective 
curvature potential is given by 



V[r(x)\ 



2M 
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In terms of the tortoise coordinate x and for x ^> M the scattering potential reads 



V(x) 



1(1 + 1) 
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The dynamics of waves propagating under the influence of logarithmic scattering poten- 
tials of the form V(x) = 1(1 + l)/x 2 + \vP xjx a (where a > 2 and f3 — 0, 1) has been studied 



in a brilliant work by Ching et al [16| 



171 ] . The study of wave tails was extended in Refs 



4C|, |4l| to include generic scattering potentials. In particular, it has been shown in 
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that the familiar case of logarithmic potentials belongs to a much wider class of scattering 
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potentials for which x^/V — > cosnt. as x — > oo. For this type of scattering potentials the 
late-time tail is dominated by 0, Q 



^~\/ {2Z) (t/2) , (5) 
where is the nth derivative of the function V(t). This simple result has one exception 
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4=11 ] : if a is an odd integer less than 2/ + 3, the late-time tail is dominated by 

v-vipm, (6) 



where V ls is the leading subdominant term in the asymptotic derivative (see |40l |41| for 
details). 

We would like to emphasize that the general expressions Eqs. (El)— (EJ) are in agreement 



with, and generalize the original results of [16|, 113] which were obtained for the specific 
case of logarithmic potentials [42]. In particular, one finds that for logarithmic scattering 
potentials the late-time tail is given by 

V^t-V^yt. (7) 

This result has the exception that if a is an odd integer less than 21 + 3, the late-time tail 
is dominated by 

* ~ r< 2,+a > In' 3 " 1 1 , (8) 

for j3 > 0, and by 

* ~ r (2/+2a-2) ^ (g) 

for = 0. 

Taking cognizance of the scattering potential Eq. (T4|), one finds that the late-time tail 
which characterizes the propagation of waves in the curved Schwarzschild spacetime is dom- 
inated by the well-known inverse power-law term 

* ~ Mr {2l+3) . (10) 
For the leading subdominant term at asymptotically late times one finds 

~M 2 t-( 2/+4 )ln 2 t , (11) 



for I > 0, and 

V l8 ~ M 2 r 4 In t, (12) 



for spherical I = perturbations. 

We would like to point out that multiple scattering from asymptotically far regions would 
also contribute to the late-time tail. For spherical perturbation (I = 0) this contribution 
decays like M 2 t~ 4 , and it is therefore negligible as compared to the leading subdominant term 
Eq. f)12p at asymptotically late-times. On the other hand, for non-spherical perturbations 
(I 7^ 0) the contribution from multiple scattering goes like M 2 H 2 ' +4 )ln 2 1. This term is of 
the same order of magnitude as Eq. (fTl~j) and it would therefore contribute to the leading 
subdominant term of the late-time tail. 

In conclusion, we have analyzed the higher-order corrections terms which "contaminate" 
the late-time tail of gravitational collapse. The existence of a relatively large correction 
term ~ M 2 t~( 2 ' +4 ) ln 2 (t/M) appears not to be widely recognized. Aside from being theo- 
retically interesting, our results are also of practical importance. They imply that an exact 
(numerical) determination of the power index (which characterizes the leading-order term 
of the late-time tail) demands extremely long integration times. (This, in turn, requires the 
use of highly stable numerical schemes.) For example, a contamination term of the form 
M 2 t~( 2l+A ^ \ia 2 (t/M) becomes ~ 1% of the leading-order term Mt~( 2i+3) at extremely late 
times of the order of 10 4 M. 
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APPENDIX A: NOTES ON A FORMER WORK 



Following the original work in Ref. |26j, Smith and Burko (SB) [36] performed an inde- 
pendent numerical computation of the late-time tail. These authors failed, however, to find 
the small logarithmic correction which accompanies the leading-order power-law tail. We 
would like to take this opportunity to make some remarks about the numerical work of SB. 



We have emphasized in Ref. [26j that it is essential to develop an extremely accu- 
rate numerical code in order to reveal the extremely small logarithmic term [which at 
late-times contaminates the (already small) leading-order power-law decaying tail]. In 



particular, in Ref. 26) we have used a modified version of the numerical scheme devel- 



oped in 
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45) . This numerical code is characterized by a fifth-order convergence 



le other hand, the numerical code 



36]. The (relatively) poor order- 



which makes it extremely stable and accurate. On t 
used by SB is only of a second-order convergence 
of-convergence of the code used by SB 36] may be the reason behind their failure to 
reveal the genuine small contamination term. 

• It should be stressed that, actually SB have reported no error in the analytical deriva- 
tion of the late-time tail presented in 26] . 



• Moreover, it is worth emphasizing that SB have not provided analytical calculations 
to support their numerical computations. (This should be contrasted with the trans- 
parent analytical derivation of the late-time tail presented in Ref. [261].) 



We would like to point out 
the results of Ching et al 



;he fact that SB provide an erroneous interpretation of 
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17| . Ching et al have shown (see also 40|, |4l|]) that for 



logarithmic scattering potentials of the form V(x) = lnx/x a with a an odd integer less 
than 21 + 3, the late-time evolution is characterized by an inverse power-law decaying 
tail without a hit factor. It should be stressed, however, that the scattering potential 



of spherical (I = 0) perturbations (the ones which were analyzed in [26l . |36| ) does 
not belong to this kind of exceptional scattering potentials. This clearly shows that 
the numerical failure of SB to detect the genuine logarithmic corrections can not be 
explained by the results of Ching et al. In fact, the results of Ching et al, both 



analytical and numerical 
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171 ]. contradict the conclusions of SB. 
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